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Abstract 

Closed form expressions for the singular-potential integrals <m\x~ a \n> are obtained 
with respect to the Gol'dman and Krivchenkov eigenfunctions for the singular potential 
Bx 2 + -^2, B > 0, A > 0. The formulas obtained are generalizations of those found 
earlier by use of the odd solutions of the Schrodinger equation with the harmonic 
oscillator potential [Aguiler a- Navarro et al, J. Math. Phys. 31, 99 (1990)]. 
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I. Introduction 

In 1990 Aguilera-Navarro et al 1 employed a perturbative scheme to provide a vari- 
ational analysis for the lowest eigenvalue of the Schrodinger operator 

d 2 A 

H= — — r + x 2 H , < X < OO, 

dx z x a 

where a is a positive constant. Writing H = Hq+XV with Hq standing for the harmonic 
oscillator Hamiltonian and V = x~ a , Aguilera-Navarro et al 1 used the basis set 

i; n (x) = A n e- x2 / 2 H 2n+1 (x), A- 2 = 2 2n+1 (2n + 1)!, n = 0, 1, 2, . . . (1) 

constructed from the normalized solutions of H ip = Eip to evaluate the matrix ele- 
ments of H. They found that 

H m+1;7l+1 = (3 + 4n)5 m , n + \<m\x~ a \n> m, n = 0, 1, 2, . . . , N - 1, (2) 

where 

^n +m >/(2n + l)!(2m + l)! 



<m|x a |n> = (— l) r 



z — n \ / 



X 

fc=0 



2 n+m n!m! 

— f. 

2 



m\ r(fc + ^)r(n - k + f ) (3) 
k) r(fc+§)r(-fc+f) 



The aim of this article is to extend these results to treat the more general spiked 
harmonic oscillator Hamiltonian 

H = --f 1 + Bx 2 + ^ + — B>0,A>0. (4) 

ax z x z x a 

The particular case of A = and B = 1 allows us, of course, to recover the results of 
Aguilera-Navarro et al 1 as a special case. The existence of all the exact eigenfunctions 
for a problem with a singular potential is the principal motivation for the present 
work: we expect that such a set of functions will be more effective for the analysis of 
other singular problems than are the Hermite functions generated, as they are, by the 
non-singular harmonic-oscillator potential x 2 . 

The article is organized as follows. In Sec. II we provide an orthonormal set of 
functions that we use to compute the matrix elements of the Hamiltonian (4) . In Sec. 
Ill we prove that this set of functions is complete in the sense that it is a basis for 
L 2 (0, oo). We then compute the matrix elements using this basis in Sec. IV. In Sec. V 
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the comparison with the result of Aguilera-Navarro et al l are presented for the special 
case B = 1 and A = 0: where we shall point out some errors in their value for the 
matrix element <ifjs\x~ a \iljs> as quoted in the appendix of Ref. [1]. 

II. An orthonormal basis 

Gol'dman and Krivchenkov 2 have provided a clear description of the exact solutions of 
the following one dimensional Schrodinger equation (in units h = 2m = 1) 



y + Vo(--- y \ ip = x E [0, oo) 

\x a J 



(5) 



with ip satisfying the Dirichlet boundary condition ip(0) = 0. They showed that the 
energy spectrum in terms of the parameters Vq and a is given by 



E n = ^VVo^n + X - + X - (Vl + 4Vbo 2 - 2a^V^j J. 



(6) 



To simplify notation we introduce the parameters B = Vqcl 2 and A = Vqo?, and obtain 
thereby an exact solution to the Schrodinger equation with the singular potential 



A 



V(x) = Bx 2 + ^r, B>0,A>0, (7) 
x z 

where the energy spectrum is now given in terms of the parameters A and B by 

E n = VB(4n + 2 + Vl + AA), n = 0,1,2,.... (8) 

The wave functions have the form 

<x\n> = ip n (x) = CnX^+^+^e-^^^i-n, 1 + -VT+4A; VBx 2 ), (9) 

2 

where n = 0, 1, 2, . . . and 1F1 is the confluent hypergeometric function 3 

The Pochhammer symbols (a)k are defined as 

(a) k = a(a + l)(a + 2) . . . (a + k - 1) = , fe = 1, 2, . . . (10) 
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where T(a) is the gamma function. Note that we have corrected the misprint of Ref. 
[4] for the power of x in the wavefunctions (9). The constant C n is determined from 
the normalization condition 

/•OO 

/ ^l(x)dx = 1, 
Jo 

which requires use of the identity 4 

oo 

n\T{v) 



^7(7+l)...(7 + n-l) 



n-l 

X 



^ ^ ^-^ n(n — 1) . . . (n — s) (7 — v — s — 1) (7 — v — s) . . . (7 — v + s) 



(11) 



[( a + l)!]2 7 ( 7 + l)...( 7 + a ) 



thereby yielding 



<2 2B5+lVi+^r(n + 1 + i VI + 4A) 



Cn n![r(l + i v / r+4Z)]2 • (12) 



In order to prove that the ip n (x), n = 0, 1, 2, ... , defined by (9) are orthonormal, we 
have to demonstrate 

00 

^ n (x)^>m(x)dx = (n^ to). (13) 
We know however, that ip n (x) and ip m (x) satisfy the Schrodinger equations 

C + (^n - Sx 2 - 4)^n = 0, £ n = ^(4n + 2 + v / TT4A), (14) 

+ (E m - Bx 2 - 4)^m = 0, E m = v / S(4m + 2 + v / l + 4A). (15) 

Multiplying (14) by ip m and (15) by ifj n and then subtracting the resulting equations, 
we obtain 

— (lpm1p' n - Ipntp'm) + ^B{n ~ m)lp n 1p m = 0. 

After integrating this equation over [0, oo) and using ip n (0) = V'm(O) = and ifj n (oo) = 
^ m (oo) = 0, we get 

oo 

(ji to) J ip n (x)ip m (x)dx = 0, 



which proves (13). Thus we obtain the following identity 
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/ 



if n ^ m, 

e- Xx \ 2 ^-\F 1 (-n, T Ax 2 )iF 1 {-m,rAx 2 )dx= { 1 n !r( 7 ) 

" if n = m, 



o l2A^( 7 ) n 
wherein the confluent hypergeo metric functions i-Fi are defined as follows 3 : 

c 

r (7) i_ 7 r£)n / 7+n _i r x = r( 7 ) 1-7 (£) _ ^n/^+n-l) /jgx 

r( 7 + n) v ; r( 7 + n) v ; v ; v ' 

for any simply closed rectifiable contour C starting at 1 and enclosing the straight 
line segment from to 1 in the complex plane, as illustrated in Fig. 1. With 7 = 
1 + \\J\ + AA and A = VB, the set of L 2 (0, oo)-functions 



1 _ 1 \„2 



il> n {x) =C n x^-^e—^ Xx iFi(-n, 7 ;A^), n = 0,l,2,... (17) 
constitutes a orthonormal system of the Hilbert space L 2 (0, oo). 

III. Proof of Completeness 

For the orthonormal functions {ip n } to qualify as a basis for L 2 (0, oo), we must 
demonstrate the density of the linear manifold generated by these functions in the 
topology induced by the norm determined by the inner product < • | ■ >. This is 
equivalent to showing that if < ip n \f >= for all n = 0, 1, 2, . . ., then / = a.e. on 
(0, oo). To this end we note that out of the fourth expression of (16) follows 

Hence the basis representation of the functions (vectors) 

UF^-n, 7 ; Ax 2 ), 1 F 1 (-(n - 1), 7; Ax 2 ), . . . , ^(-1, 7; Ax 2 ), ^(-0, 7; Ax 2 )} 
in terms of the basis 

{x 2 ™,^- 1 ),...,^ 2 ,!} 
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is achieved by a lower triangular (n + 1) x (n + 1) matrix, whose diagonal entries are 
r(7+n-fc) f° r ^ = 0, 1, 2, . . . , n - i.e. this matrix is invertible provided A ^ 0. Thus 
each x 2n is a unique linear combination of the n + 1 functions i-Fi(— (n — k), 7; Ax 2 ) for 
fe = 0, 1, 2, . . . , n, which conclusion carries over to the 2n-th degree Taylor polynomial 



fix 2 ^ 1 / ^r 2 \ 
fc=0 v 7 



of e 2 about the point 0, where \i is an arbitrary parameter. 

Let / be an L 2 (0, oo)-function orthogonal to each of the ip n , which is equivalent to 
saying 



00 

< e n (-^ )\f>=J x^-h n (-^f)f(x)dx = (19) 





for all n = 0, 1, 2, . . .. Here we note that 



x 1 2 e 4 j( x ) 

in an L 1 (0, oo)-function whose absolute value majorizes x 1 ~^e~ h ^ e^^ f{x) for < 
\ and consequently also x^~^e~ A ^e n (— fix 2 / 4) f (x) . 

Because jT"3e"T"e„(- fix 2 /4)f(x) converges to x 1 ~^e~ 2 ^~f(x) a.e. on (0, 00) as 
n — > 00, we conclude by means of the Lebesgue dominated convergence theorem 5 that 
we may replace e n (— [xx 2 /2) by in Eq.(19) for all complex numbers \x such that 

I A* I < 1? which, after setting x = y/2t, yields the Laplace- Transform expression 

£{F}(z) = J e~ zt (V2iy-%f{V2i)dt = 0, \z - A| < j. (20) 


However, the Laplace transform of the measurable Laplace-transformable function 
F(t) = e _zt (v / 2t) 7_ ^/(v / 2t) defines a holomorphic function of variable z in the right 
half plane $t(z) > vanishing in the disc \z — A| < ^. By uniqueness of the analytic 
function 6 the Laplace transform of the function must vanish in the right half plane, 
specifically £{F}(s) = for all s on the interval (0, 00). Further, the Laplace trans- 
form determines F(t) uniquely 7 a.e. in t on (0, 00), hence F(t) = a.e. in t or / is the 
zero L 2 (0, oo)-function. Consequently, {-i/Vi : n = 0, 1, 2, . . .} is an orthonormal basis of 
L 2 (0,oo). 
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IV. The matrix elements <m\x a \n> 

Let us now split the Hamiltonian (5) into an Hq part 



d 2 A 
ff = - +fe 2 + a;>0 (21) 



dx 2 x 
and a perturbation 

Hi = \/x a . 

The eigenstates of Hq are now given by (9) and their unperturbed energy is given by 

(8) . All we need to do is to evaluate the matrix elements < m\x~ a \n > using the basis 

(9) , namely 



< mix a 



oo 

\n >=C n C m J e-^ x2 x-« +1+VT + TX 1 F 1 (-n, 1 + ^VT+4A; VBx 2 ) 



x iFi(-m, 1 + -Vl + 4A; VBx 2 )dx, 
2 



(22) 



a < 2 + y/l + 4A 
This is equivalent to 

< m\x~ a \n >= 9j^ B -^- a+2+ ^A) x 7j (23) 

where 

oo 

7 = y e~ r r 7_s iFi(-n,7; r)iFi(-m, ^;r)dr (24) 
o 

with r = v 7 ^ 2 , 7 = 1 + ±Vl + 4A, and s = 1 + § . 

^From the Fubini-Tonneli theorem 5 combined with the Leibniz formula for dif- 
ferentiating the product of two functions, as well as exponential shift from the third 
expression to the fourth in Eq.(16), with n replaced by m, we find that / is given by 
the expression 

7-n-l / -(l-t)r l-s 



V ; r(n + 7 )r(m + 7 ) V 1 J y } J 

c o 

^(-l) fc ( ™ ) ( 7 + m - 1)( 7 + m - 2) . . . (7 + fc) r 7+m-i-(m-fc) 
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Further, owing to the fact that the simply closed rectifiable contour C lies to the left 
of the complex number 1, Fig. 1, we have 

oo 

J e -(l-t)r r J-8+k dr = T ^_ s + k + 1 ^ 1 _ t yy+s-k-l 


and our expression for / thereby reduces to 

T _ ( n n n![r( 7 )] 2 " fc /m\ r(m + 7 )r( 7 -s + fc + i) 
1 } r(n + 7 )r(m + 7 )^ J U; r(fc + 7 ) 

k ~° (25) 



x (2ni)- 1 ^t- n - 1 (l-t) 

c 



Since the contour C" has in its inside, Fig. 1, and the integrand has a weak singularity 
at 1 (in consequence of 3?(s + n — k — 2) > —1), the Cauchy integral formula lets us 
write the contour integral multiplied by (27u) -1 as the n-th derivatives of (i_^ s + n - fc - 2 
evaluated at t = 0. Utilizing thereafter the Pochhammer symbol in Gamma function 
format Eq.(10), we arrive at 



T Bf fv ,/m\ r( 7 - g + Hl)r( g + n-fc-l) 

' r(n + 7 )^ 1 ' \k) r( 7 + jfe)r(s - k - 1) 



Therefore, the matrix elements are given by 

™ fc /ro\ r(fc + l + Iy^+4A-f)r(f +n- 

W r(fc+i + WT+4Z)r(f -a>) 



2 



X 

/— 

a < 2 + Vl + 44, 



with normalization coefficients C n given in Eq.(12). In case ^ — /c is a negative integer, 
then 3 l/r(^ — /c) = for such /c and the terms involving these fc's shall not appear 
in the summation of Eq.(27). Further, by expressing the confluent hypergeometric 
functions i-Fi(— n, 7; r) and i-Fi(— m, 7 ;r) by means of the fourth formula in Eq.(19) 
and substituting these into Eq.(18), we immediately see that the sum appearing in 
Eq.(27) is a polynomial of degree m + n in a. 
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With Eq.(27) we have therefore computed the matrix elements of the operator 
x~ a in the complete basis given by the Gol'dman and Krivchenkov eigenfunctions (9). 
Concomitant to our result, the matrix elements 



<U\x- a \n> = {-l) n B^ 



T{l + yi + AA) r(l + Wl + 4A-f)r(f + 



n 



n\T(n + l + yi + 4A) T(l + + 4A)r(§ ) 

(28) 

are of special interest. 

V. Explicit forms of the matrix elements 

In terms of the parameter 7 = 1 + ^y/T+ 4~A, the explicit forms of the first ten matrix 
elements of x~ a are for a < 27: 

a _„ f rX-f+7) 
x oo —° r( 7 ) 

a _ gf IH+7) 

x oi - 2 Vrr( 7 ) a 

_bt r(-f+ 7 ) , , 

x ° 2 ~ 22 ^^Tr )r{l) a ^ + 2 > 

03 23 v /3!7(7+l)(7+2)r( 7 ) V 7 

*n =S-S^(a 2 -2a + 47) 

a"" =-fff r (-f+^) «(«2 _ 2a + g ) 

12 8 v ^+I)r( 7 +i) v ,; 

x7 3 a =Sf . r ^ +7 > a( tt 3 + (127 - 4)a + 24 7 ) 

13 i6 v /3!( 7 +2)( 7 +i)r( 7 +i) v v / y 

=^ f W^g} (« 4 " 4« 3 - 16(1 + 2 7 )a + 4(3 + 4 7 )a 2 + 32 7 (1 + 7)) 

X ™ 32V3i 2 K7+ 2 )r(7+2) a(a4 " + 4(5 + 67) " 2 " 16(2 + 37)tt + %7(1 + 7)) 

X " a 3 r 8 ( 4 ~^ 7 3 ) ) (a 6 - 8a 5 + (72 + 367) a 4 - (208 + 144 7 )a 3 + (272 + 72O7 + 

2887 2 )a 2 - (192 + 1152 7 + 5767 2 )a + 384 7 (1 + 7) (2 + 7)). 
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These matrix elements can be compared for the special case (A, B) = (0, 1) with the 
matrix elements computed by the simple harmonic oscillator representation supple- 
mented by Dirichlet boundary condition [1]. We have found an error in the value of 
the matrix element as given by [1]; this error is confirmed by a re-computation 
according to the matrix element expression given by them. Indeed, the matrix element 
£33" should read: 

x -£ = y r V (« 6 - 6a 5 + 106a 4 - 384a 3 + 2080a 2 - 3408a + 5040) 
7.1 (2) 

instead of 

= _j_2Z( a 6 _ 6a 5 + 106a 4 _ 454a 3 + 1660a 2 _ gggg^ + 504() ) 
7.1 



as quoted by Aguilera-Navarro et al 1 . 

In this work we have proved that the Gol'dman-Krivchenkov wavefunctions con- 
stitute an orthonormal basis for the Hilbert space L 2 (0, 00). Using this orthonormal 
basis, we are able to construct the matrix elements of x~ a : the general result (Eq.27) 
is convenient for use in any practical application which involves such singular potential 
terms. It is also interesting that, with minor changes 8 , essentially involving only the 
value of the coefficient A, the same formulas apply immediately to the corresponding 
problems with non-zero angular momentum and in arbitrary spatial dimension N > 2. 
A detailed variational analysis of the spiked harmonic Hamiltonian operator based on 
these matrix elements is presently in progress. 
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Figure 1 The Contour C" in the t-plane starting at 1 and enclosing the line segment 
from to 1. 



